Adiabatic Magnetization of Superconductors as a High-Performance Cooling 

Mechanism 
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The adiabatic magnetization of a superconductor is a cooling principle proposed in the 30s, 
which has never been exploited up to now. Here we present a detailed dynamic description of the 
effect, computing the achievable final temperatures as well as the process timescales for different 
superconductors in various regimes. We show that, although in the experimental conditions 
explored so far the method is in fact inefficient, a suitable choice of initial temperatures and metals 
can lead to unexpectedly large cooling effect, even in the presence of dissipative phenomena. Our 
results suggest that this principle can be re-envisaged today as a performing refrigeration method 
to access the fiK regime in nanodevices. 
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Since the very early discovery of the laws of ther- 
modynamics, cooling represents one of the most fasci- 
nating challenges for both experimental and theoretical 
physics [H, 0). A well known cryogenic principle is the 
adiabatic demagnetization, based on the property that 
ordinary magnetic materials, such as paramagnetic salts, 
experience an entropy decrease when a magnetic field is 
applied, due to the alignment of their atomic dipoles. 
This effect is currently applied also to nuclear spins un- 
der large magnetic fields, allowing to reach the /xK regime 
in nuclear demagnetization refrigerators 

In superconducting materials, however, the opposite 
cooling principle is observed. It is well known [3[ that 
a sufficiently strong magnetic field drives a superconduc- 
tor (S) into the normal (N) state, and that such phase 
transition occurs with a supply of latent heat, since the 
S state is much more ordered phase than N at a given 
temperature. As a consequence, if a magnetic field with 
intensity H increasing from up to the critical value H c 
is quasi-statically applied on a thermally isolated super- 
conductor, its entropy S per unit volume is preserved 



S N (T f ,H = H c )^S s (T t ,H = 0), 



(1) 



and the metal cools from the initial temperature Tj 
down to a final temperature Tf [1, 0], as illustrated in 
Fig. [TJ This cryogenic principle, known after the pio- 
neering works by Mendelssohn & Moore Q and by Kee- 
som & Kok ,7] in 1934 as "adiabatic magnetization of 
a superconductor" (AMS), offers the advantage that the 
required magnetic fields are much lower (H <C IT) than 
those typically used in adiabatic demagnetization refrig- 
erators. Furthermore, the presence of a metal greatly 
simplifies the contacting to devices and allows faster equi- 
libration time-scales. Although the validity of AMS was 
successively confirmed by other experiments [1, 0, Hcj |. 
only a relatively small cooling effect was observed so far: 
the temperature was lowered from 2.50 K to 2.22 K on 




FIG. 1: (color online) Adiabatic magnetization of a super- 
conductor. Solid curves describe the entropy S of a metal in 
the N and in the S state as a function of the reduced tem- 
perature. When a magnetic field is applied on a thermally- 
isolated superconductor at T%, the metal is driven into the N 
state {A — * B), and the temperature decreases down to Tf. 
Dashed lines refer to the entropy in the intermediate state. 



tin samples [a, [llj, from 1.43 K to 1.32 K on thallium 
samples, and from 3.63 K to 3.54 K on lead spheres [Io| . 
It thus never became of practical use as a cryogenic tech- 
nique, and its theoretical modeling has also been over- 
looked. On the other hand, the exponential growth of 
nanotechnological applications at low temperature de- 
mands higher performance to refrigerators, which are re- 
quired to be more versatile, faster and not invasive. The 
aforementioned features of AMS seem quite promising to 
this purpose, and a detailed analysis of this refrigeration 
principle is desirable. 

Here we present the first dynamical description of the 
adiabatic magnetization effect, taking into account the 
role of dissipative phenomena, and computing both the 
final temperature and the process time-scales. This anal- 
ysis allows us to show that, while the conditions of the 
experiments carried out so far were not suitable for cool- 
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ing, realistic regimes can be identified in which the adi- 
abatic magnetization can be exploited as a performing 
cooling technique. 

We consider the case of type-I superconductors, and 
start our analysis with some remarks about thermody- 
namics. In each phase the entropy includes a phonon 
and an electronic contribution, S N ^ S = S p h + S^/ S . Ex- 
plicitly, S p h(T) = aT 3 , where T is the temperature and a 
is the coefficient related to the Debye temperature. The 
electronic entropy in the N state has a linear behavior 
S^(T) — jT. The contribution of spin paramagnetism 
is negligible in the range of magnetic fields we are inter- 
ested in, H < IT, so that <S N (T, H c ) S 5 N (T, 0). In the 
S phase the condensate is a coherent state with vanishing 
entropy, so that 5| ; is purely due to quasi-particles and 
can be obtained from the BCS theory as 
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FIG. 2: (color online) Final temperature Tf vs initial temper- 
ature Ti for several type-I superconductors, in the absence of 
dissipative effects. 



S S el (T) 



dEN(E,T) [f(E) ln(f(E))] , (2) 



where vp is the normal density of states (DOS) 
at the Fermi level, f(E) = (1 + exp[^/(fc B T)])- 1 
is the Fermi-Dirac distribution function, Af(E, T) = 
\E\j^jE 2 - A(T) 2 Q(E 2 - A(T) 2 ) is the BCS normalized 
DOS, with A(T) denoting the superconducting order pa- 
rameter and O(x) the Heaviside function. 

For a given initial temperature Tj, the final tempera- 
ture Tf of the metal is determined by Eq. |T]) , which can 
be rewritten as 



T, 



f 



Tf 



(T*) 2 (T*) 



T 3 
_ % 

*\2 



(3) 



indicating that Tf depends in general on two character- 
istic parameters, namely the critical temperature T c of 
the superconductor, and T* = 



h/a = ^hZTy^ 2 T Fl 
which defines the temperature below which the entropy 
of the N state is dominated by the electron contribu- 
tion. Here, Z denotes the nominal valence, while Tp 
and Td the Fermi and Debye temperatures of the metal, 
respectively. Furthermore, $ is a universal function 
of T/T c defined through the relation S% (T)/S ph {T) = 
(T*/T c ) 2 $(T/T C ), exponentially small for T/T c < 0.1, 
and of the order of unity for 0.5 < T/T c < 1. Despite the 
simplicity of its derivation, Eq. ([3]) contains important 
physical insight. Indeed if the initial temperature T is 
of the same order as T* , so is the final temperature Tf 
( T f ^Ti), even if T < T c . In this regime the AMS is 
therefore clearly inefficient as a cooling mechanism. By 
contrast, if T <C T*, the final temperature decreases as 



Tf ~Tf/(T*f 



(4) 



This cubic dependence stems from the fact that, in this 
temperature regime, the AMS effectively transforms the 
entropy of phonons into the entropy of electrons. Wc 
emphasize that this effect represents an advantage with 



respect to the linear gain factor Tf /T characterizing the 
adiabatic demagnetization process The efficiency of 
the AMS thus heavily depends on the material choice and 
on the initial temperature range. 

Figure [2] shows Tf vs T calculated for several super- 
conductors from 10 mK to the zero-field critical temper- 
ature T c . The low-temperature linear behavior in the 
log-log plot accounts for the cubic dependence ([!]), and 
occurs in all materials. At higher temperatures, differ- 
ences emerge between metals that exhibit T* < T c (like 
Pb, Hg, Tl), and those with T* > T c (like Al, Zn, and 
Ta). In the former case Tf < T, whereas in the lat- 
ter case Tf exhibits a steep decrease governed by T c , i.e. 
Tf ~ Tf <5>(T/T C )/T 2 , before the crossover to the cubic 
law ([4]). It is noteworthy that most experiments were con- 
cerned with the first group of metals, and with ~ T*. 
This explains the unsatisfactory cooling observed on Sn 
d, [n[ and Pb [l(| ■ Our analysis suggests that tantalum 
(Ta) is a good candidate, since T*/T c > 3, and it allows 
to obtain Tf in the range ~ /iK...mK starting from T in 
the range ~ 100mK...l K. Notice that Sn is suitable only 
if T < 0.6 K, whereas Al is even better for T < 0.2 K. 

So far, using purely thermodynamical arguments 
and BCS theory, we have shown that AMS may in 
principle lead to extremely low values of Tf, provided 
the superconducting metal and initial temperature 
are appropriately chosen. However, the influence of 
dissipative effects must be taken into account in order for 
such a method to be considered as a promising cooling 
technique. First of all, when a magnetic field is applied 
to a type-I superconductor, the transition to the N phase 
is preceded by formation of an intermediate state (IS), 
where S and N phases coexist for H'JT) < H < H C (T). 
Here H' C (T) = (1 — n)H c (T), n is the demagnetization 
factor, H C (T) is the critical field defined through the 
relation dH 2 c {T)/dT = {2/^ )[S s (T, 0) - S N (T, 0)] 0], 
and /iq is the vacuum permeability. The presence of 
a normal fraction a; at in the IS yields dissipative eddy 
currents when the magnetic field is increased with 
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time. Furthermore, in a cryostat the superconductor is 
connected to some mounting support that remains at 
the initial temperature X^; the metal is thus exposed 
to a heat flux, which cannot be neglected in view 
of its relatively small low-temperature specific heat. 
Finally, once the cooling is realized, heating generated 
by measurements on any device attached to the cryostat 
has to be considered. Thus, even assuming that the 
range of initial temperatures and the superconductor 
are properly chosen, the existence of dissipative effects 
leads to the following questions: i) is the AMS-based 
cooling robust against these effects? ii) if so, what are 
the typical time-scales in which low temperatures are 
reached, and how long can these be maintained? 

To address these questions, we have analyzed the AMS 
dynamically, i.e., the process is governed by the equation 
d t S = P(t)/T(t), where P is the total dissipated power 
per unit volume which involves the three contributions 
mentioned above. For simplicity, we consider a bundle 
of N w long and thin superconducting wires of radius R 
and length L each, attached to an insulating support of 
length I, as sketched in the inset of Fig. 3(a). The AMS is 
driven by the magnetic field, and three time regimes can 
be distinguished. In the first one the magnetic field is in- 
creased from to H' c (Ti), and no cooling occurs since the 
whole system remains superconducting, so that P = 
and T = Tj. In the second regime (cool-down) H is var- 
ied from H' c (Ti) up to H c (0) over a time r, and the system 
enters into the IS state. Although a detailed description 
of the IS for a given geometry is, in general, quite compli- 
cated, we wish here to capture its main characteristics. 
We shall thus follow Ref. [5j by assuming that the N and 
S regions of the IS are uniformly distributed, so that the 
normal fraction xn increases with the magnetic field as 
x N (T,H) = 1 - n _1 (1 - H/H C (T)), and that the en- 
tropy is a linear combination of the N and S entropies, 
S(T,H) = x N S N {T,0) +(l-x N )S s (T,0). In this case 
the AMS is described by the differential equation 

Cv (T 1 H)T~^T ( ^±H = P , (5) 
n dl 

where 

C v (T,H)^x N C^(T)+(l-x N )C^(T)+C^(T,H) (6) 

is the total specific heat (per unit volume) in 
the IS, Cy /S (T) = TdS s / N {T)/dT, and C l f = 
(TH/Li nH*(T))(S N (T, 0) - S S {T, 0)) 2 is the latent heat 
to be supplied for the transition. As soon as xm ^ the 
cooling mechanism is enabled, and the temperature of the 
metal starts to lower, though contrasted by the dissipa- 
tive effects. A simple calculation shows that the variation 
of the magnetic field B = xn[1qH c (T) trapped in the nor- 
mal fraction induces eddy current dissipation P e ddy(t) = 
R 2 <tB 2 /8 in each wire, where a is the electric conduc- 
tivity of the metal. At the same time, the temperature 
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FIG. 3: (color online) (a) Time evolution of Tj in the in- 
termediate state calculated at Ti = 250 mK for three values 
of t. The inset shows a scheme of the adiabatic magnetiza- 
tion cooler, (b) Full time evolution of Tj calculated at three 
different T for t = 1 h. Dashed curves refer to Pi oa d 7^ 0: 
from bottom to top, P io(ld = 10 pW, 1 nW, and 100 nW. All 
calculations were performed for Ta. (see text) 



gradient across the insulating support between the 'hot' 
upper surface at temperature Ti and the 'cold' lower sur- 
face in contact with the metal [see the inset of Fig. 3(a)] 
leads to a heat flow P supp = fe(lf +1 - T p+1 )/((3 + 1)IL. 
Here, b and (3 are parameters characterizing the tempera- 
ture dependence of thermal conductivity n supp (T) = bT@ 
of the insulating support 1] . Any temperature gradient 
in the metal has been neglected due to its relatively high 
thermal conductivity. 

Figure 3(a) displays the temperature evolution for 
< t < t calculated for three values of t, starting from 
Ti = 0.25 K. For simplicity we assumed H to vary lin- 
early with time. For t <C t the temperature experiences 
a relative slow decrease, whereas a fast drop is observed 
for t < r. The cooling effect is eventually contrasted by 
both dissipative eddy currents and heat flow from the 
support. It is worth emphasizing that P e ddy and P S upp 
behave differently with respect to the velocity of the mag- 
netic field variation. For fast field variations [solid curve 
of Fig. 3(a)] Peddy is relevant and P SU pp is suppressed. 
In contrast, for slow field variations [dotted curve of Fig. 
3(a)] Joule heating has a minor effect while the heat flow 
from support affects the cooling for longer time. For a 
given geometry and initial temperature, the competition 
of these two terms determines the optimal time which 
allows to reach the lowest temperature, as shown by the 
dashed curve in Fig. 3(a) for r = 1 h. Such time-scale de- 
pends on the electric conductivity of the metal, and the 
thermal conductivity of the support. Tantalum seems to 
be a good candidate superconductor due to its relatively 
low conductivity a ~ 10 9 £l~ x vnT x , and high specific heat 
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(7 ~ 523Jm- 3 K- 2 and a ~ 2.63 Jm- 3 K- 4 ). We note 
that aluminum (Al), in spite of its high ratio T*/T c , is 
less suitable for AMS due to its extremely high electric 
conductivity; alternatively, tin (Sn) may be a fair choice. 
As far as the support is concerned, a good insulator like 
PVC, with parameters b = 1.8 X 1CT 5 Wm^K" 1 and 
f3 = 2.05 seems appropriate. The plots of Fig. 3 refer 
to this case 

The last time-regime corresponds to the case where 
the system is fully normal. The magnetic field is not fur- 
ther varied [H = H c (0)], so that P e dd y = 0. This is the 
regime where measurements are typically carried out on 
a device thermally anchored to the metal. We have thus 
included a constant load power for t > t arising from 
the measurement (Pi oa d) besides P sup p. The result of the 
whole dynamical process is shown in Fig. [3fb) for three 
initial temperatures T, corresponding to the base tem- 
perature of a 4 He cryostat (Ti = 1.2 K), a 3 He cryostat 
(Ti = 250 mK), and a dilution refrigerator (T { = 20 mK). 
For each T, the solid (dashed) curve represents the tem- 
perature evolution without (with) Pi oa d- Notably, this 
cooling method ensures in all these ranges of operation 
a temperature gain of about two orders of magnitude, 
which can be reached within an hour or less, and can be 
maintained for several hours. This represents an advan- 
tage with respect to the time-scales typical of the adia- 
batic demagnetization of nuclei. The external load sus- 
tained by the AMS method depends on the temperature 
range of operation. For instance, we have calculated that 
a bundle of wires of about 4 x 10 3 cm 3 volume operating 
at T — 1.2 K can sustain a power load of 100 nW with- 
out significantly affecting its final temperature, whereas 
at T — 20 mK a power load of 10 pW increases Tt of 
few hundreds of after 5 hours of operation [see Fig. 
3(b)]. We stress that ordinary superconducting electron- 
ics (such as tunnel junctions circuits, radiation detectors, 
and SQUIDs) as well as single-electron devices exhibit 
power dissipation typically below 1 pW, thus suggesting 
that AMS is suitable to operate on nanostructures in the 
ultra-low temperature regime. 

Finally we notice that, since the magnetization must 
evolve through equilibrium states, the variation of the 
applied magnetic field must proceed slow enough for re- 
laxation processes to ensure equilibrium between elec- 
trons and lattice phonons. The determination of such 
time-scales in the IS is a crucial issue, since the N and S 
phases have much different characteristic relaxation rates 
[14j . Analyzing the three terms of Eq. ([B]) , one can easily 
prove that even a small normal fraction xn ~ 10~ 3 is 
sufficient for Cy to largely dominate the other two con- 
tributions. Thus, apart from an extremely small range of 
magnetic fields, the specific heat of the superconductor in 
the IS is essentially determined by the electronic contri- 
bution in the N fraction, which drives the cooling " drag- 
ging" the lattice phonons and the S fraction. An upper 
bound for the relaxation time characterizing the process 



is therefore represented by the inverse of the electron- 
phonon scattering rate in the N phase, which scales as 
T el-ph x t 2], and is typically much shorter than that 
of the superconducting phase. For Ta in the temper- 
ature range 2 x 10 -4 . . . 5 x 10 _1 K, T e i- p h lies in the 
. 10 3 : 



14| , thus ensuring the consistency 



range ~ 10 
of our quasi-static approach. 

In conclusion, we have provided a dynamic descrip- 
tion of cooling by adiabatic magnetization of supercon- 
ductors. We have shown that, while in the experimental 
conditions explored so far the method is in fact inefficient, 
a suitable choice of temperature ranges and superconduc- 
tors make this principle promising as a high-performance 
refrigeration technique. Beside involving low magnetic 
fields (i.e., ~ 10~ 2 to 10 _1 T), the present method of- 
fers the additional advantage that the final temperature 
depends cubically on the initial one [see Eq.j4])]. More- 
over, we find that the cool-down times are comparable or 
shorter than those of typical demagnetization cryostats in 
the same temperature range, while the warming-up rates 
can be of the order of several hours under continuous 
power load (see Fig. [3]). Our results suggest that magne- 
tization cycles to improve this cooling principle can also 
be envisioned. 
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